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Transport through two one-dimensional interacting met- 
als (Luttinger liquids) coupled together at a single point is 
analyzed. The dominant coupling mechanism is shown to 
be of electrostatic nature. Describing the voltage sources by 
boundary conditions then allows for the full solution of the 
transport problem. For weak Coulomb interactions, transport 
is unperturbed by the coupling. In contrast, for strong inter- 
actions, unusual nonlinear conductance laws characteristic for 
the correlated system can be observed. 

PACS numbers: 71.10.Pm, 72.10.-d, 73.40. Gk 

The physics of one-dimensional (ID) conductors has 
received much attention lately, chiefly due to fabrication 
advances and the discovery of novel ID materials such as 
carbon nanotubes From the theoretical point of view, 
these systems are of interest since Coulomb interactions 
invalidate the ubiquitous Fermi liquid description. The 
resulting state is often of Luttinger liquid (LL) [||J|] type 
characterized by, e.g., spin-charge separation, suppres- 
sion of the tunneling density of states, and interaction- 
dependent power laws in the transport behavior. How- 
ever, so far the unambiguous experimental observation of 
LL behavior has been difficult to achieve. 

In this paper, we study two correlated ID metals cou- 
pled in a point-like manner ("crossed Luttinger liquids"). 
For the standard two-chain problem, where two Luttinger 
liquids are connected all along the conductors, the cou- 
pling normally destroys the LL phase [Q. In the case 
of a point-like coupling, however, the LL characteristics 
can survive and lead to the unusual transport features 
reported below. The most promising candidates for their 
experimental observation are carbon nanotubes. At not 
exceedingly low temperatures, metallic single-wall nan- 
otubes exhibit LL behavior (with an additional flavor 
index) ||. In a remarkable recent experiment, Tans et 
al. Q were able to attach leads to a single nanotube. 
So far, transport measurements have been dominated by 
Coulomb charging effects due to rather large contact re- 
sistances between the leads and the nanotube, thereby 
masking any possible deviation from Fermi liquid the- 
ory. In the near future this problem might be overcome, 
and non-Fermi liquid laws should indeed emerge. Other 
realizations of crossed Luttinger liquids could be based 
on, e.g., ID quantum wires in semiconductor heterostruc- 
tures , or edge states in a fractional quantum Hall bar 




FIG. 1. Two Luttinger liquids coupled together at one 
point (x = 0) and connected to external reservoirs held at 
constant voltages ±f/i/2 and ±(72/2. In the absence of sin- 
gle-particle tunneling, the currents obey 7^ = I[. 

The geometry of our system is shown in Fig. |^. We 
shall consider two spinless Luttinger liquids characterized 
by the same interaction constant g |^ . Here the noninter- 
acting value is 5 = 1, and externally screened Coulomb 
interactions imply < 5 < 1 For the nanotube 

experiment of Ref. |^ , one has an externally unscreened 

/\x — x'\ interaction potential and therefore very strong 
correlations. Strictly speaking, this interaction leads to 
(7 in an infinite system, but the finite length of the 
nanotube in Ref. |^ implies g « 0.2. A natural and 
quite simple description of Luttinger liquids is offered by 
the standard bosonization method js). External reser- 
voirs (voltage sources) can be incorporated by imposing 
boundary conditions |10[| for the phase fields employed in 
the bosonization scheme. This approach offers a general 
and powerful route to studying multi-terminal Landauer- 
Biittiker geometries for strongly correlated electrons. 
The crossed Luttinger Hquids depicted in Fig. |] may be 
the simplest example for such a problem. 

We start by expressing the right- and left-moving 
{p = ±) component of the electron operator ippi{x) in 
conductor « = 1 or 2 in terms of the dual bosonic phase 
fields Oi{x) and (j)i{x) obeying the algebra 

[Mx),e,{y)]- ^ -{t/2)6,,sgn{x - y) . (1) 

The bosonization formula then reads (3| 

ippi{x) = IL- cxpl-ipkpx ~ ips/TTg9i{x) - i^n / g (f)i{x)] 
V 27ra 

(2) 
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where the same average density ^f/tt is assumed for both 
conductors. The short-distance cutoff (lattice spacing) in 
Eq. (^) is taken as a = l/fej?- To ensure anticommutation 
relations among different branches {p,i), we use (real) 
Majorana fermions ripi fulfilling [ripi,ripiii]+ = 2Spp'Siii. 
In the following, only products of Majorana fermions will 
appear. A valid choice for these products employs the 
standard Pauli matrices ||^, 



VplV-p,2 = 



(3) 



Assuming that the conductors do not contain impurities, 
the Hamiltonian of the uncoupled system is 



(4) 



where we have put Ti — \ and the sound velocity v = 
vp/g = 1. Adiabatically connected voltage sources can 
then be taken into account by Sommerfeld-like boundary 
conditions. Applying the voltage Ui along conductor 1, 
and U2 along conductor 2, see Fig. |], they read ||l0|] 



{Pp=±,i{^ too)) = ± 



eUj 
Airg 



(5) 



where p± i is the density of right/left-moving particles 
injected into conductor i. Outgoing particles are assumed 
to enter the reservoirs without reflection. 

Let us now consider a point-like coupling of both ID 
conductors at, say, a: = 0. For example, in a nanotube 
setup, two nanotubes could be stacked on top of each 
other. Such a contact causes (at least) two different cou- 
pling mechanisms JT^ . 

First, there arises a (density-density) electrostatic in- 
teraction of the form Vi ~ Aipi(0)/02(0). Using Eqs. (||) 
and (^, and omitting the mean density /cf/tt which is 
supposedly neutralized by positive background charges, 
the bosonized representation of the density operator is 



i{x) =p — sm[2kpx - 
TTa 



^0,{x)] , (6) 



where the first ( "slow" ) term is due to the sum of right- 
and left-moving densities pm + pLi, and the second 
("fast") term arises from mixing right- and left-movers. 
The =F signs correspond to i ~ 1,2, respectively. One 
checks easily that most contributions to Vi are irrele- 
vant for g < 1, i.e., they have scaling dimension 77 > 1. 
Keeping the fast component in Eq. (^ yields the only 
important term. 



Vi = - 



Ai 



(tta)' 



■sin[v/4^^i(O)]sin[v/4^02(O)] , (7) 



with scaling dimension 771 = 2g. Clearly, this cou- 
pling becomes relevant for sufficiently strong interactions. 



g < 1/2. In contrast, a static potential scatterer in one 
of the conductors would be relevant already for 5 < 1 
]l3[ . In our case, electrons in conductor 1 experience 
the fluctuating potential scattering Vi due to electrons in 
conductor 2, implying a doubled scaling dimension. 

The second potentially important process is single- 
particle hopping from one conductor into the other. It is 
helpful to distinguish processes that do (do not) preserve 
the p = R, L ^ ± index, yielding the two perturbations 
V2 = A2 J2p V'pi(0)V'p2(0)+ H.c. (preserving the p index) 

and V3 = A3 Y,p V'pi(O)'0-p2(O)+ H.c. (not preserving the 
p index). They have the bosonized form 



2A 

V2 = -<J, cos{^ [01 (0) - 62(0)]} 



xsin{V^[0i(O)-02(O)]} 



2A 



V3 = —<JySin{^[9i{0) + 92m} 
na 



(8) 
(9) 



For the standard two-chain problem, a (bulk) coupling 
term formally identical to Eq. (H) has been discussed 
in Ref. iQ. Both V2 and V3 have scaling dimension 
V2 = Vi ^ (5 + 1/5)72 > Ij from which one might 
naively conclude that they are irrelevant [ p^ . However, 
this conclusion is premature because V2 and V3 have 
conformal spin 5* = 1 |16|. For an operator with non- 
zero conformal spin, the standard criterion for relevance 
T] < 1 does not apply, since relevant perturbations may be 
generated in higher orders of the renormalization group 
(RG). This phenomenon indeed occurs in the standard 
two-chain problem [ pT[, where the (bulk) coupling term 
corresponding to Eq. (|8|) generates relevant particle-hole 
and/or particle-particle excitation operators. Similarly, 
we find that V2 and V3 together generate the electro- 
static coupling Vi given in Eq. (Q), but no other relevant 
terms. Omitting irrelevant operators, the resulting RG 
equations take the closed form 



^ = [l-2g]\+[g-l/g]t\ 
| = [l-(5 + l/5)/2]^ 



(10) 



where A2 = A3 = i is the hopping amplitude and Ai = A 
the electrostatic coupling. The standard flow parameter 
is defined by di = —dlnuc, where uJc is a high-frequency 
cutoff that is reduced under the RG transformation. 

Let us first discuss the case g = 1. The electrostatic 
coupling is irrelevant, i.e., we may effectively put A = 0, 
but the hopping term stays marginal. By refermioniz- 
ing the Hamiltonian Hq + V2 + V3, and employing the 
boundary conditions (^, one arrives at the familiar re- 
sults for uncorrelated electrons in the geometry of Fig. 
see Ref. Q- We therefore recover the usual Landauer- 
Biittiker formalism. Second, for g < 1, the hopping am- 
plitude t{£) always scales to zero as ^ ^ 00, and the 
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effects of single-particle tunneling can be captured by a 
renormalization of the bare electrostatic coupling A, see 
Eq. (0). We shall assume henceforth that this renormal- 
ization has been carried out, and only the electrostatic 
interaction Vi will be kept. In that case, the currents 
flowing through conductor i = 1 or 2 satisfy — I-, see 
Fig. |l|, and can be computed from the bosonized current 
operator |H| 



= ey 



0,i) 



(11) 



For weak interactions, 1/2 < g < 1, the electrostatic 
coupling X{£) also flows to zero as ^ ^ oo. In that case, 
at low energy scales, crossed Luttinger liquids are ba- 
sically insensitive to the coupling considered here. At 
asymptotically low energy scales, the currents are then 
li — {e^/h)Ui. The finite-temperature or low-voltage 
corrections due to the irrelevant operators Vi can be com- 
puted by perturbation theory in the respective coupling 
strengths A^. Since the fluctuating potential scattering is 
irrelevant for 1/2 < 51 < 1, the corrections due to V2.3 are 
governed by the standard exponent a = {g -\- g~^ — 2)/4 
for tunneling into a bulk LL (HI). This is in contrast to 
a static potential scatterer, where tunneling into the end 
of a LL matters at low energy scales . 

Directly at g = 1/2, the operator Vi is marginal, and 
straightforward refcrmionization yields 



h 1 + (A/27ra): 



(12) 



Each conductor exhibits a response only to the voltage 
applied to itself, with the conductance now explicitly de- 
pending on the electrostatic coupling strength A. 

For sufficiently strong interactions, g < 1/2, the elec- 
trostatic coupling A(£) flows to strong coupling. To pro- 
ceed, we switch to the linear combinations 



d±{x)^{ei{x)±92{x)}/V2, 
^±(x) ={(/)i(a;)±(/)2(a;)}/V2 , 



(13) 



which again obey the algebra (|l|). Remarkably, the 
Hamiltonian Hq + Vi decouples into the sum H+ + H- 
with 



H± = \j dx{{d^^±f + {d,d±f] 
^gd±{Q) 



(14) 



± 



■ cos 



and the boundary conditions (g) determining the density 
pp^r of p = ± movers injected into channel r = ± take 
the form 



{Pp,r{x -poo)) 



p e{Ui+rU2) 
V2 47r5 



(15) 



Therefore we are left with two completely decoupled 
systems r = ±, each of which is formally identical to 



the problem of an elastic potential scatterer embedded 
into a spinless LL ||l^. However, this LL now has the 
doubled interaction strength parameter g — 2g. The 
boundary conditions ( [l5| ) specify the effective voltages 
Ur = {Ui + rU2)/ \/2 applied to channel r = ±. In anal- 
ogy to Eq. (p^ ) , currents in channel r = ± are defined by 
Ir = e^Jg/n dt'dr{0), and from Eq. (p^, we then find the 
currents h = (J+ ± /_)/ V2 flowing in conductor i = 1,2. 

The Hamiltonian dl4) has been discussed in detail be- 
fore, see, e.g., Refs. p^l8|-^. For arbitrary g, the exact 
solution of the transport problem has been given in Ref. 
p9[ . This solution exploits the integrability of Eq. ( p^ ) 
and employs the thermodynamic Bethe ansatz. Simpler 
exact solutions are possible by means of refermionization 
techniques for 5 = 1 [see Eq. (p^] and g = 1/2. The case 
(7 = 1/2 thus corresponds to an uncorrelated situation in 
the new basis (|l3|), and g = 1/4 is the Toulouse point 
pO| . Progress can also be made by expanding in |e| <C 1 
forg=l-e ]l| org = l/2-e p. 

Employing the exact results of, e.g., Ref. at zero 
temperature we find the asymptotic low- voltage behavior 

- X— I sgn(C/i + U2) [e\Ui + U2\/\b]'/'-' (16) 

± sgn(C/i -U2)[e\Ui- U2\/ Xb]'^'-'} , 

where the ± sign corresponds to i — 1,2, respectively. 
The energy scale As generated by the bare electrostatic 
coupling A is given by 



Xb = {cg/a) (A/a)i/[i-29] , 



(17) 



where Cg is a numerical constant of order unity p9[ | . The 
result (n6h holds under the condition 



e\Ui±U2\ « Xb ■ 



(18) 



If both voltages approach zero, the linear conduc- 
tance vanishes in both ID conductors. We thus find 
a pronounced zero-bias anomaly, with characteristic 
interaction-dependent power laws for small voltages. 

Let us now discuss the full current-voltage character- 
istics. A particularly simple solution emerges at the 
Toulouse point g — 1/4 by refermionization ||2^ of 
Eq. ( p^ ) under the boundary conditions (15). At zero 
temperature, the result is 



J± = {e^h) [U± - V±] 



(19) 



where V± is the four-terminal voltage [|lO| subject to the 
self-consistency equation 



eV± = 2ABtan-i{[2e;7± - i3/2)eV±]/ Xb} 



(20) 



where Xb — A^/4(7ra)'^ in accordance with Eq. (p^. Un- 
der the condition (18), the exact result ( p9[ ) reproduces 
Eq. ( p^ ) again. In the absence of a coupling, Xb = 0, one 



finds the correct unperturbed currents li 



yh)u,. 
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FIG. 2. Transport current Ji normalized to the unper- 
turbed value /? = {e^/h)Ui for g = l/4,As/e = l,r = 0, 
and several values of the cross voltage U2- 

The transport current Ji is plotted as a function of 
Ui in Fig. ||. Contrary to what is found in the uncor- 
related system fll|| , the current Ii is extremely sensitive 
to the applied cross voltage U2- For U2 — 0, transport 
becomes fully suppressed for Ui —^ 0, with (7-dependent 
nonlinear low- voltage corrections given by Eq. (16). In- 
creasing U2 for some fixed Ui then leads to an increase 
in the current Ii . Eventually, the linear conductance be- 
havior is restored for very large cross voltage. In fact, for 
e\Ui±U2\ ^ As, one always recovers the unperturbed 
currents Ii — {e'^/h)Ui. The generic correlation effects 
are most important under the conditions (p^. 

Remarkably, there is a suppression of the current if 
Ui = ±1/2, which is observed as a "dip" in the nor- 
malized current displayed in Fig. ^ This effect can be 
rationalized in terms of a partial dynamical pinning of 
charge density waves in conductor 1 due to commensu- 
rate charge density waves in conductor 2. As can be 
checked from Eq. ([l9|), while the nonlinear conductance 
Gil = dli/dUi stays positive, one can have a nega- 
tive value for 6*12 = dIi/dU2- In fact, the latter (off- 
diagonal) conductance is within the bounds —e^/2h < 
G12 < e^/2h [we note that G12 = for g > 1/2], while 
the diagonal conductance fulfills < Gn < e'^/h. The 
pronounced and nonlinear sensitivity of the current Ii 
to the applied cross voltage U2 is a distinct fingerprint 
for Luttinger liquid behavior. Parenthetically, anoma- 
lous power laws can also be found in the temperature 
dependence of the current. 

To conclude, we have examined nonlinear transport 
through two Luttinger liquids coupled together at one 
point. The only relevant coupling is of electrostatic ori- 
gin and leads to distinct correlation effects for strong 
Coulomb interactions, g < 1/2. The theoretical find- 
ings reported here could be of use for the experimen- 
tal identification of non-Fermi-liquid behavior in carbon 
nanotubes and other one-dimensional materials. 

We thank A.O. Gogolin, H. Grabert, and C.A. Stafford 
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